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ABSTRACT 

Context. We study the sharpness of the time-resolved prompt emission spectra of gamma-ray bursts (GRBs) observed by the Gamma- 
ray Burst Monitor (GBM) on board the Fermi Gamma-ray Space Telescope. 

Aims. We aim to obtain a measure of the curvature of time-resolved spectra that can be compared directly to theory. This tests the 
ability of models such as synchrotron emission to explain the peaks or breaks of GBM prompt emission spectra. 

Methods. We take the burst sample from the official Fermi GBM GRB time-resolved spectral catalog. We re-fit all spectra with a 
measured peak or break energy in the catalog best-fit models in various energy ranges, which cover the curvature around the spectral 
peak or break, resulting in a total of 1,113 spectra being analyzed. We compute the sharpness angles under the peak or break of the 
triangle constructed under the model fit curves and compare them to the values obtained from various representative emission models: 
blackbody, single-electron synchrotron, synchrotron emission from a Maxwellian or power-law electron distribution. 

Results. We find that 35% of the time-resolved spectra are inconsistent with the single-electron synchrotron function, and 91% are 
inconsistent with the Maxwellian synchrotron function. The single temperature, single emission time, and location blackbody function 
is found to be sharper than all the spectra. No general evolutionary trend of the sharpness angle is observed, neither per burst nor for 
the whole population. It is found that the limiting case, a single temperature Maxwellian synchrotron function, can only contribute up 
to 58;^jg% of the peak flux. 

Conclusions. Our results show that even the sharpest but non-realistic case, the single-electron synchrotron function, cannot explain 
a large fraction of the observed GRB prompt spectra. Because any combination of physically possible synchrotron spectra added 
together will always further broaden the spectrum, emission mechanisms other than optically thin synchrotron radiation are likely 
required in a full explanation of the spectral peaks or breaks of the GRB prompt emission phase. 

Key words, gamma rays: stars - (stars:) gamma-ray burst: general - radiation mechanisms: non-thermal - radiation mechanisms: 
thermal - methods: data analysis 


1. Introduction 

The prompt emission of gamma-ray bursts (GRBs) is one of the 
most puzzling observed astronomical phenomena. Since the dis¬ 
covery of GRBs in 1967 (Klebesadel et al. 1973), many emis¬ 
sion models have been proposed in order to explain the prompt 
phase of gamma-ray emission. This phase consists of gamma 
rays mainly within tens to a few hundred keV, in some cases as 
high as a few thousand keV, lasting from a few milliseconds to 
hundreds of seconds. 

Gamma-ray bursts are distributed isotropically (Briggs et al. 
1996; Hakkila et al. 1994; Tegmark et al. 1996) and cosmolog- 
ically (Metzger et al. 1997; Waxman 1997) over the sky. De¬ 
spite the last 45 years of research efforts, the dominant emission 
mechanism of these cosmological sources is still controversial. 
Synchrotron radiation from a simple electron population is one 
of the simplest physical phenomena that may be able to produce 
the observed spectral slopes of the Band function (Band et al. 
1993) that is commonly used to describe the photon spectra of 
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GRB prompt emission. The Band function is an empirical math¬ 
ematical function consisting of two segments of power laws, de¬ 
scribed by the low- and high-energy photon indices a and fi, con¬ 
nected at the peak energy parameterized as Ep. This peak energy 
has been observed typically at hundreds of keV (Kaneko et al. 
2006; Nava et al. 2011; Goldstein et al. 2012, 2013; Gruber et al. 
2014). In what is known as the fireball model (Goodman 1986; 
Meszaros et al. 1993; Meszaros & Rees 1993; Rees & Meszaros 
1992, 1994; Tavani 1996; Piran 1999), there are ejected shells 
with different bulk Lorentz factors. When the faster shells catch 
up with the slower shells, internal shock waves will be produced. 
The electrons in the shocked region of the shells are accelerated 
and their energy is radiated via synchrotron emission in the local 
magnetic field. 

The Band function’s two power-law indices are usually com¬ 
pared to the slopes of various radiation models, leading to the 
discovery of the so-called line-of-death problem (Katz 1994; 
Crider et al. 1998; Preece et al. 1998, 2002; Tavani 1995) for the 
synchrotron theory. When a power-law distribution of electron 
energies is combined with synchrotron radiation theory, the low- 
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energy power-law photon index is -2/3. The fact that a fraction 
of observed a > -2/3 indicates that, at least in some cases, the 
synchrotron explanation of GRB prompt spectra can be problem¬ 
atic, because the observed spectra rise faster. The observed vio¬ 
lations of the line-of-death are typically around 30%. Recently, 
Burgess et al. (2015) used the Band function to fit a large num¬ 
ber of simulated slow-cooling synchrotron spectra (with spectral 
peak determined by injection energy of the electrons rather than 
their energy losses, as in the fast-cooling case), concluding that, 
in practice, the line-of-death may be steeper, a > -0.8, than the 
value of a > -2/3. Moreover, they found that the Band function 
cannot recover the simulated synchrotron peaks and power-law 
indices. These findings all question the validity of using the syn¬ 
chrotron theories to explain the Band parameters. 

Instead of fitting the empirical Band function to the spec¬ 
tra, Burgess et al. (2014) used a synchrotron function in the fit¬ 
ting process, combining the slow-cooling scenario with thermal 
emission. Yu et al. (2015) used a double broken power law to 
fit eight bright GRBs, in which they found that the line-of-death 
problem could be alleviated in a moderately fast-cooling sce¬ 
nario, in which the fast- and slow-cooling electrons are mixed 
together, usually with a blackbody component at tens of keV 
or in a varying magnetic field. However, no single synchrotron 
model could completely explain all the spectral properties of 
GRB prompt emission. 

In this work we study the sharpness of the synchrotron emis¬ 
sion spectrum in comparison to time-resolved spectra of GRBs, 
a question recently raised by Beloborodov (2013); Vurm & Be¬ 
loborodov (2015). Our approach focuses on the curvature of the 
region capturing the peak or break energy in the GRB prompt 
spectra by re-fitting all the spectra in an energy domain depend¬ 
ing on this peak or break, using the burst sample from the Fermi 
GBM time-resolved spectral catalog (Yu et al. in prep.). By com¬ 
paring the spectral sharpness of the observed spectra to vari¬ 
ous physical emission models, we are able to directly determine 
whether a model is capable of accounting for the peak emission 
of the observed spectra. By concentrating on the spectral peak or 
break, we avoid potential issues with interpretation of the asymp¬ 
totes of the fit functions, which might lie outside the observable 
domain or be contaminated by additional radiative processes or 
instrumental effects. 

Thanks to the high-quality gamma-ray data obtained by 
Fermi GBM (Meegan et al. 2009), which provides wide energy 
coverage and fine temporal and spectral resolutions, this is the 
first time that we can directly compare the curvature of a large 
number of time-resolved spectra to that of physical models, so 
that statistically significant conclusions about the prompt emis¬ 
sion mechanism can be drawn. 

This paper is structured as follows. In Section 2, we describe 
our analysis method. The results are presented in Section 3. In 
Section 4, we check the consistency of our analysis. In Section 5, 
we discuss the theoretical implications. The summary and con¬ 
clusions are given in Section 6. Unless otherwise stated, all er¬ 
rors reported in this paper are given at the Icr confidence level. 

We note that a recent independent study of the peak-flux 
GRB spectra (Axelsson & Borgonovo 2015) shows that a syn¬ 
chrotron function could be too wide for the observed Band 
shape. They measured the full-width-half-maximum in the vFy 
spectra obtained from the 4-years Fermi GBM GRB time- 
integrated spectral catalog (Gruber et al. 2014; von Kienlin et al. 
2014) and the BATSE 5B GRB spectral catalog (Goldstein et al. 
2013). 


2. Data analysis 

2.1. The data and the method 

The Fermi GBM consists of 12 thallium activated sodium iodide 
(Nal(Tl)) detectors, which cover 8 keV - 1 MeV, and 2 bismuth 
germanate (BGO) detectors, which cover 200 keV - 40 MeV. The 
combined energy range of the two kinds of detectors is ideal for 
the study of GRB prompt emission spectra because the typical 
spectrum peaks at a few hundred keV. 

In order to account for the change in orientation of the source 
with respect to the detectors, due to the slew of the spacecraft, 
RSP2 files are used in the fitting process, which contain the de¬ 
tector response matrices for every two degrees on the sky. For 
each burst a low-order polynomial (order 2 - 4) is fit to every 
energy channel, according to a user-defined background inter¬ 
val before and after the prompt emission phase, and interpo¬ 
lated across the emission interval. All spectra are re-fit with the 
GBM official spectral analysis software RMFIT' v4.4.2BA and 
the GBM response matrices v2.0. 

Our sample is taken from the official Fermi GBM GRB time- 
resolved spectral catalog (Yu et al. in prep.) which consists of the 
brightest bursts observed by GBM before 21 August 2012. All 
the bursts in our sample are long bursts, i.e., with Tgo > 2 s (Kou- 
veliotou et al. 1993). They were selected according to 3 criteria: 

(1) the total fluence in 10 keV - 1 MeV, / > 4.0 x 10“^ erg cm“^; 

(2) the peak flux in 10 keV - 1 MeV, F^ > 20 ph s“' cm“^ in 
either 64, 256, or 1,024 ms binning timescales; and (3) the burst 
has 5 or more time bins when binned with signal-to-noise ratio 
SIN - 30. This results in 81 bursts and 1,802 spectra in total, 
of which 311 do not satisfy the catalog’s standard error criteria 
(for the details on the error criteria, see Gruber et al. 2014). We 
exclude from further analysis these 311 spectra, and concentrate 
on the remaining 1,491 spectra. 

Only the spectra best fit by the Comptonized model (COMP), 
the Band function (BAND), and the smoothly broken power law 
(SBPF) are included in the analysis. The functional forms of 
these models are given in Appendix A for completeness. This 
is because we are interested in comparing the sharpness around 
the peak or break energies of theoretical models to the observed 
spectra. Thus, the 194 spectra best fit by a simple power law or 
the power law plus blackbody, are excluded. All the 1,297 spec¬ 
tra best fit by either COMP, BAND, or SBPF have convex shape 
(i.e., a > —2 for COMP and a > j3 for BAND and SBPF). 

The best-fit model parameters for the 1,297 spectra are ob¬ 
tained from the catalog. Using the catalog values of Fp, every 
spectrum is re-fit (using the same best fit function as in the cat¬ 
alog) in a narrower energy domain that covers Ep (or the break 
energy Eb if there is no peak in the spectrum). We refer this en¬ 
ergy domain as the "data domain", which contains the "triangle 
domain" (described below). We find that in the 1,297 spectra, 34 
of them have no converged re-fit and are therefore excluded. In 
the remaining 1,263 spectra, 150 of them have large error bars 
(according to the criteria from Gruber et al. 2014) and thus are 
further excluded. In total, 1,113 spectra are used in this work, of 
which 942 are best fit by COMP, 99 by BAND, and 72 by SBPF. 

Motivated by the necessity to model the spectral curvature 
around Ep and exclude any possible curvature contribution from 
the low- or high-energy tail, we test the data domain on a few of 
the brightest bursts (details are discussed in Sect. 4.2). As a re¬ 
sult, (Fieft, Fright) = (O.lFp, 3.0Fp) is adopted for a good balance 
between statistics and optimal description of the spectral sharp- 

' The public version of the RMFIT software is available at http: // 
fermi.gsfc.nasa.gov/ssc/data/analysis/rmfit/ 
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Fig. 1. Illustration of how the triangle is constructed and the sharpness 
angle 6 is defined. The shaded regions indicate the data domain (see 
Sect. 4.2). The triangle is constructed within the triangle domain (see 
Sect. 4.3), under the best-fit model curve (black). The vertical and hor¬ 
izontal axis are plotted in logarithmic scale in units of normalized vFy 
flux and photon energy, respectively. 


ness. Assuming a typical spectram with Ep ~ 300 keV, it means 
that we are covering the range from ~ 30 keV to ~ 900 keV. 

The peak energy Ep and the peak flux vFy(Ep) are used to 
normalize the model curve such that the peak coincides with 
{x,y) = {EIEp,vFy{E)lvFy(Ep)) - (1,1). For the spectra of 
SBPL fits without a peak, the break energy E^, is used instead. 
For each spectrum, a triangle ((1,1), (xmi, t/ieft), (-bright, l/right)} be¬ 
low the spectral curve is constructed in dimensionless space as 


f (1,1), ] 

= \ (£'ieft/£'p,vFy(£'ieft)/vFy(£'p)), i. (1) 

[ (£'i-ight/£'p,vFv(£'iight)/vFv(£'p)) j 

The sharpness angle 6 is computed in logarithmic space, under 
(log 1, log 1) and between log xieft and log xught (which we refer 
to as "triangle domain"). Thus, 6 is an indication of the spectral 
sharpness and is independent of the actual position of Ep (i.e., 
also independent of redshift). Figure 1 illustrates how such a tri¬ 
angle can be constructed. 

Similarly, we also construct the right-angled triangle 
{(1,1), (xieft, 1/left), (1,1/left)) below the spectral curve between xieu 
and X - 1, and compute the left-hand side angle flieu in logarith¬ 
mic space, under (log l,log 1) and between logxieft and log 1. 
Thus, in the limit of small xieu, 0ieft becomes equivalent to a mea¬ 
sure of the steepness of the low-energy power-law slope. 


2.2. Synchrotron emission modeis 

The monochromatic flux of the synchrotron emission spectrum 
can be obtained, given the electron population n^, as 


Fy OC 



«e(re) T 



dye. 


where 


T{x) = X 



( 2 ) 


(3) 


in which is the Lorentz factor of the electron, Ve is the syn¬ 
chrotron frequency of the electron, and Ks/^ is the modified 
Bessel function of fractional order 5/3. For a single electron, 
the synchrotron spectrum is simply proportional to 'F.^ It can 
be shown that the limits of F(x) can be approximated by sim¬ 
pler analytical functions for x « 1 and x » 1 (see, e.g., van 
Eerten & Wijers 2009), for the ease of computation. Notice that 
Eqn. (2) either describes an instantaneously generated spectrum 
and 90 degrees pitch angle between magnetic held and electron 
velocity, or a situation where magnetic held and particle popula¬ 
tion remain unchanged. 

Mathematically, the synchrotron emission spectrum of a sin¬ 
gle electron is the sharpest case. However, under realistic condi¬ 
tions, there is no reason to believe that the observed data orig¬ 
inates from only one electron. Thus, it is more realistic to con¬ 
sider a Maxwellian population of electrons, since it is an efficient 
distribution of electron energies and sharper than typical non- 
thermal spectra. Eor a Maxwellian population of electrons with 
the temperature parameterized by the thermal Lorentz factor yth, 
we have 


rie 



(4) 


and the Maxwellian synchrotron spectrum 



Since Ve oc y^, by changing the variable v = ^Vth oc it can 
be shown that 


J X 5exp(-^5;c 2)T(x)dx, (6) 

Tih 

which allows us to normalize the spectrum in units of Again, 
we note that Eqn. (5) represents one of the sharpest cases among 
synchrotron spectra for multiple electrons, but that the assump¬ 
tions of a single temperature and magnetic held are still unreal¬ 
istic. Observed emission will contain a mixture of these and lead 
to smoother spectra. 

Another reasonable assumption for the electron population 
is a power-law distribution of the electron energies; 

Ue oc yl’’ : ye > ym , (7) 


and the synchrotron spectrum with population index p is 

^00 

Fycx y-/r 
dym 

where y^ is the minimum injection energy of the electron popu¬ 
lation. The Fy spectrum can be solved as 




(9) 


where is the minimum injection frequency of the electron 
population. As with temperature, the observed spectrum will be 
smoother due to a mixture of values in the emission. If we 

^ This already assumes integration over emission direction (see Ry- 
bicki & Lightman 1986, Eqns. 6.29 - 31). If a single electron were 
viewed from a single angle, a sharper spectrum would mathematically 
result. 
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Fig. 2. Left panels: Cumulative distribution functions of 6 and distributions of <rg. Right panels: Cumulative distribution functions of Sieft and distri¬ 
butions of CTieft. The limits of the normalized blackbody (dotted line), single-electron synchrotron (solid line), and synchrotron with a Maxwellian 
distrihution function (dashed line) are overlaid. In the above legends, COMP represents the Comptonized model, BAND represents the Band 
function, SBPL represents the smoothly broken power law, and ALL represents the overall population (COMP + BAND -l- SBPL). 


substitute the approximation of 'Fix) ~ for x <s: 1, we can 
recover the 1/3 low-energy slope below Vm for any value of p. 

In reality, electron cooling should exist, as more energetic 
electrons lose their energy faster due to radiative losses and cool 
down. One could consider, in addition to the minimum injec¬ 
tion energy break Vm, the cooling break Vc (see, e.g.. Fig. 1 in 
Yu et al. 2015). However, the ratio between Vm and Vc depends 
sensitively on assumptions on the shock micro-physics and fluid 
evolution. Additionally, the sharpness of the cooling break de¬ 
pends intrinsically on the distribution of electrons throughout the 
shock region, and has no local analog for a single electron pop¬ 
ulation. For the purpose of this paper, it is sufficient to consider 
the case without cooling: quick evolution of electron energies 
due to cooling will smoothen the synchrotron spectrum. There¬ 
fore, any cooling synchrotron spectrum can never be sharper 
than Eqn. (8). 

For comparison, we also consider blackbody emission, 
which is given by 


oc 


exp(hvlkT) - 1 


( 10 ) 


3. Spectral sharpness results 

Figure 2 (left panels) shows the cumulative distribution functions 
(CDFs) of the sharpness angles 6 and the distributions of the er¬ 


Table 1. Sharpness angle 9 and left angle 6ieft for various emission mod¬ 
els. “We note that if p < 3, vFy keeps on increasing monotonically. 


Emission models 

6 (degrees) 

0ieft (degrees) 

Blackbody 

43 

27 

Single-electron synchrotron 

97 

53 

Maxwellian synchrotron 

135 

64 

Synchrotron with p = 2“ 

170 

40 

Synchrotron with p = 4 

128 

56 


rors cTg. The dotted, solid, and dashed black vertical lines indi¬ 
cate the values of 9 for the normalized blackbody, single-electron 
synchrotron emission function,^ and synchrotron emission func¬ 
tion from a Maxwellian electron distribution, from left to right. 
These values are listed in Table 1. It is found that over 35% of the 
spectra are inconsistent with single-electron synchrotron emis¬ 
sion and 91% are inconsistent with synchrotron emission from 
a Maxwellian electron distribution. The blackbody spectrum is 
found to be much sharper than any of the observed spectra. 

^ For a single emission direction, 9 and flieft are 76 and 43 degrees re¬ 
spectively, for the polarization direction perpendicular to the projection 
of the magnetic field on the sky, and are 67 and 37 degrees respectively, 
in the parallel case. These values reflect the standard textbook results 
for single-electron emission prior to convolving with an electron distri¬ 
bution function (see Rybicki & Lightman 1986). 
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The synchrotron emission function from a Maxwellian elec¬ 
tron distribution produces one of the sharpest (i.e., narrowest) 
spectra (Sect. 2.2). The values of 6 for the synchrotron emission 
function from a power-law electron distribution for p - 2 and 
p = 4 are also listed in Table 1. The spectrum for p - 2 was 
normalized by the peak position in the Fy space, because for 
p < 3, vFy keeps on increasing monotonically. Notice also that 
the spectrum for /? = 4 is of similar sharpness to the spectrum 
for Maxwellian (for p ^ oo, the spectrum would reduce to a 
single-electron synchrotron spectrum). 

In principle, erg should be propagated directly from the errors 
on the observed photon counts, since the counts are independent 
of the choice of the fitting models. However, the spectral peak 
can only be found and the flux can only be normalized when the 
counts are convolved with a model (COMP, BAND, etc.) and 
the response matrices, through RMFIT. Therefore, we compute 
CTg by performing Monte-Carlo simulations using the errors of 
the re-fit model parameters. First, we extract the Icr etTors from 
the RMFIT results. Because the errors on model fit parameters 
a and p (see Appendix A) are not necessarily Gaussian, we then 
randomly draw new values of a and p from a uniform probability 
function sharing the same Icr error, and we re-compute 9. This 
process is repeated 1,000 times for every spectrum. We then take 
the Icr width of the resulting 9 distribution and average over left 
and right Icr values. We note that our method generates the most 
conservative values of erg, since the uniform probability function 
has the largest standard deviation. 

As shown in the bottom left panel of Fig. 2, the resulting dis¬ 
tribution of CTg has a median around 5 degrees. This is too small 
to affect our conclusions. However, we note that erg for BAND 
and SBPL can be systematically larger than those for COMP, be¬ 
cause the high-energy tail of COMP is an exponential cutoff with 
no parameter dependence. Therefore, a has very little effect on 
the right-hand-side spectrum for COMP fits (see Eqn. A.l), and 
CTg of COMP may be under-estimated. 

Because the fit results for BAND and SBPL fits are dis¬ 
tributed over a wide range and have larger angles and errors than 
COMP, it is of interest to look separately at the low-energy left- 
hand-side angles flieft. This way we can explore how both sides 
contribute to the total curvature and shape our results. Also, 0ieft 
is unaffected by the transition from photon counts to upper lim¬ 
its that sometimes already occurs slightly below S.OLp on the 
right-hand-side. We therefore show, in the right panels of Fig. 2, 
the CDFs of the low-energy left-hand-side angles 6>ieft (i.e., the 
angle under (log 1, log 1) and between log ;icieft and log 1) and the 
distributions of their errors erjeft. 

The top right panel of Fig. 2 shows that if one compares 0ieft 
instead of 9, the overall fraction inconsistent with single-electron 
synchrotron increases to 48%, and the overall fraction inconsis¬ 
tent with Maxwellian synchrotron is also 91%. Therefore, even 
when only the left-hand-side of the spectral peak (or break) is 
considered, the same conclusions can be drawn, and the errors 
remain sufficiently small not to affect the final result. The dis¬ 
tributions of the errors on 0ieft are similar for the different fit 
functions, suggesting that their values are not merely driven by 
the curvature of the fit function itself. In addition, it shows that 
the low-energy curvature is the main cause of the violation of 
any synchrotron emission model, and that the upper limits in the 
high-energy side could harden the high-energy power laws of 
BAND and SBPL, which make the spectral shape less sharp. 

Of our 1,113 spectra, 35% violate the synchrotron line-of- 
death (i.e., a > -2/3), higher than the 20% observed by Gruber 
et al. (2014) in their peak-flux "P" spectra sample. This implies 
that a large number of spectra are still consistent with the line- 



0 20 40 60 80 100 

maximum fraction of Maxwellian at x = 1 (%) 

Fig. 3. Distribution of the maximum fraction contributed from the 
Maxwellian synchrotron function at x = 1. The solid histograms rep¬ 
resent the distrihutions using the best-fit model parameters, while the 
dashed histogram shows the minimum allowed sharpness by the uncer¬ 
tainties from the best-fit parameters. Spectra with 100% at x = 1 are 
accumulated in the last bin. 



Fig. 4. Example spectrum taken from GRB 101014.175 (2.560 - 
3.584 s), showing the maximum contrihution to the best-fit model 
by the Maxwellian synchrotron function, at x = 1. The normalized 
Maxwellian synchrotron (green curve) and the best-fit model (black 
curve) overlaid. The black dashed lines show the peak position of the 
best fit model and the relative normalized flux levels. In this particular 
spectrum, the Maxwellian fraction is about 65% at x = 1. The shaded 
regions show the boundaries Xieft = 0.1 and = 3.0. Deep green 
data points are from the BGO detector and the others are from the Nal 
detectors. Triangles represent upper limits. For display purpose, the bin 
size has been increased by a factor of 5 - 10 relative to the standard hin 
size. 


of-death. However, we find that in the 65% of spectra that do not 
violate the line-of-death, 92% of them violate the Maxwellian 
limit (i.e., 0 < 135 degrees) given in this paper. This shows that 
the sharpness angle method can identify many more spectra that 
are consistent with the line-of-death but are still sharper than 
what the synchrotron theory predicts. By contrast, of the 35% of 
spectra that violate the line-of-death, only 10% of them do not 
violate the Maxwellian limit. 
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Fig. 5. Six examples of evolutionary trends of 6. Red, blue, or green color indicates that the best-fit model is COMP, BAND, or SBPL, respectively. 
The light curves are overlaid in arbitrary units. The limits of the normalized blackbody (dotted line), single-electron synchrotron (solid line), and 
synchrotron emission from a Maxwellian electron distribution (dashed line) are overlaid. 


Since our results indicate that the synchrotron model alone 
cannot explain most of the prompt spectra, we can ask the ques¬ 
tion, if synchrotron emission is still one of the mechanisms that 
contributes to the observed peak flux, how much at most can it 
realistically contribute? In Fig. 3, we show the distribution of 
the maximum peak flux contributed from the Maxwellian syn¬ 
chrotron function. For the spectra that do not have a peak, we 
compute this value at the spectral break. A sample spectrum 
from GRB 101014.175"^ is plotted in Fig. 4. The normalized 
Maxwellian synchrotron function was shifted vertically and hor¬ 
izontally until the distance between its value at x = 1 and the 
peak of the fit model is minimized. The advantage of evaluating 
this value at the peak of the fit model is that it is energy domain 
independent. For the spectra of SBPL without a peak, the break 
energies E\, are used instead. It is found that the Maxwellian can 
only contribute up to 58;^jg% of the peak flux (solid histogram). 
Even if the minimum sharpnesses (i.e., the broadest) allowed by 
the uncertainties in the best-fit parameters are considered, this 
only slightly increases to 68 ^ 23 *^ (dashed histogram). Again, we 
caution that these synchrotron spectra represent a limiting case 
of high sharpness, relative to that expected from a distribution of 
temperatures and magnetic held strengths and a rapidly evolving 
particle population. In that sense, 58% indicates an upper limit. 

3.1. Spectral evolution 

We now consider the sequence of spectra within bursts. We se¬ 
lect and plot in Fig. 5 the evolution of 0 for 6 example bursts, 
with the Maxwellian synchrotron limit and the observed light 
curves overlaid. It can be seen that 8 exhibits various evolution¬ 
ary trends; 

In this paper, the names of the bursts are given according to the Fermi 
GBM trigger designation that is assigned for each new trigger detected. 
The first 6 digits indicate the year, month, and day of the month, and 
the last 3 digits indicate the fraction of the day. For more details, please 
see the online Fermi GBM burst catalog at http://heasarc.gsfc. 
nasa.gov/W3Browse/fermi/fermighrst.html 


(1) In GRB 100414.097 (top left panel), the spectrum be¬ 
comes less sharp as time evolves. We also plot the spectra of this 
burst in Fig. 6 , with the normalized blackbody (red), Maxwellian 
synchrotron (green), and the best-fit model (black) overlaid. 
The violation of the Maxwellian synchrotron function is clearly 
shown in this example, and 9 increases with time. In the typical 
fireball model, 6 is expected to increase with time due to, e.g., 
increasing collision radius and curvature effects. 

(2) In GRB 110731.465 (top central panel), the opposite hap¬ 
pens and the spectrum becomes sharper as time evolves. 

(3) In GRB 120711.115 (top right panel), 9 fluctuates be¬ 
tween the limits of single electron and Maxwellian, without clear 
correlation to the observed light curve. We note that the first time 
bin at around the trigger time has a small 9. 

(4) In GRB 090902.462 (bottom left panel), 9 remains ap¬ 
proximately constant in the plateau during the first 7 s, and then 
increases to higher but fluctuating values (11 - 25 s). We note that 
during 7 - 11 s, the catalog best-fit model is the power law plus 
blackbody, in accordance with the finding of Abdo et al. (2009). 
We did not compute the sharpness angle for this period of time 
because the blackbody is sharper than all synchrotron cases. 

(5) In GRB 090926.181 (bottom central panel), the low emis¬ 
sion level first time bin gives the largest 9, which is consistent 
with what Maxwellian synchrotron emission predicts (in con¬ 
trast to GRB 090902.462), and 9 then decreases and fluctuates 
around the value of the single-electron limit. It increases again 
in the penultimate time bin to a value marginally consistent with 
the Maxwellian limit, and then drops again to the single-electron 
limit. 

( 6 ) GRB 090829.672 (bottom right panel) has the largest 
fraction of spectra consistent with a Maxwellian synchrotron 
explanation (13 out of 32 spectra, 40%). Similar to GRB 
090926.181, it combines a large value of 9 with a low emission 
level. During the main emission pulses between 35 - 55 s, 0 de¬ 
creases below the Maxwellian limit and then increases again to 
values above the Maxwellian limit. 
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Fig. 6. Spectral evolution of GRB 100414.097 with the normalized blackbody (red), Maxwellian synchrotron (green), and the best-fit model 
(black) overlaid. Time evolves from top left to bottom right, and the time since trigger is labeled at the top of each snapshot spectrum, in units of 
seconds. The peaks of the models are all normalized to {x,y) = (1,1). Data points and the shaded regions are plotted as described in Fig. 4. For 
display purpose, the bin size has been increased by a factor of 5 relative to the standard bin size. 
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Fig. 7. Left panel: Cumulative distribution functions of 9 for the catalog full energy domain fits. Right panel: Distributions of ag. The limits of the 
normalized blackbody (dotted line), single-electron synchrotron (solid line), and synchrotron emission from a Maxwellian electron distribution 
(dashed line) are overlaid. In the above legends, COMP represents the Comptonized model, BAND represents the Band function, SBPL represents 
the smoothly broken power law, and ALL represents the overall population (COMP -h BAND + SBPL). The blue dotted line and histogram show 
fit results if all spectra are fit using BAND, provided that they are converged fits, but not necessarily the best fit when compared to other models. 


These bursts are chosen to show the variety of evolutionary 
trends in 0: gradual increase, gradual decrease, fluctuation be¬ 
tween the single-electron and Maxwellian limits, small 6 during 
low emission level and large 0 during high emission level, large 
9 during low emission level and small 0 during high emission 
level, and decrease from above the Maxwellian limit followed 
by an increase again to above the Maxwellian limit. 


4. Consistency checks 

4.1. Choices of the fitting modeis 

It is observed that over 66% of the time-resolved catalog best-fit 
models are COMP. This indicates that most of the observed spec¬ 
tra are indeed sharper than BAND or SBPL would predict. The 
same statistical behavior was also observed in the GBM time- 
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Fig. 8. Comparison of the convolved data points and the respective 
convolving model curves for a sample spectrum taken from GRB 
100414.097. The red curve and data points are obtained from the COMP 
fit, the blue ones are from the BAND fit, and the orange ones are from 
the SBPL fit with the break scale A allowed to vary. The Maxwellian 
synchrotron function is also overlaid (green). For display purpose, the 
bin size has been increased by a factor of 5 relative to the standard bin 
size. 


integrated spectral catalogs (Goldstein et al. 2012; Gruber et al. 
2014) and the BATSE time-integrated spectral catalogs (Kaneko 
et al. 2006; Goldstein et al. 2013). 

We show in Fig. 7 the CDFs of 6 and the distributions of erg 
from the catalog best fits, evaluated using the full data domain of 
(8 keV, 40 MeV) and triangle domain of (xieft, xnght) = (0.1,3.0). 
It is observed that the catalog best fits produce results similar to 
the re-fits. 

We note that COMP is inherently an exponential cutoff 
model, while BAND and SBPF are power laws joined by a peak 
or break energy. This intrinsic difference between the fit func¬ 
tions motivates us to explore the fit results if all spectra are fit 
using the Band function. Therefore, we further plot in Fig. 7 
the distributions of 0 and erg using the catalog BAND fits ("all 
BAND"), provided that it is a converged fit with a peak (or break) 
energy in the vFy space, but not necessarily the best fit when 
compared to other models. We find that it gives larger 9 and erg. 
This indicates that when the Band function is applied to all spec¬ 
tra, the values of 9 can be over-estimated due to larger uncer¬ 
tainties. Nevertheless, even in the all-BAND approach, 77% of 
spectra are sharper than the Maxwellian synchrotron limit. 

The models of COMP, BAND, and SBPF have been ex¬ 
tensively tested over the years and are found to provide good 
fits to data (e.g., Kaneko et al. 2006). In Fig. 8, we show the 
comparison of the convolved data points and the respective con¬ 
volving model curves for an illustrative sample spectrum taken 
from GRB 100414.097 (see also Fig. 6). The red curve and 
data points are obtained from the COMP fit (CSTAT/dof^ = 
301.66/285), the blue ones are from the BAND fit (CSTAT/dof 
= 301.61/284), and the orange ones are from the SBPF fit 
(CSTAT/dof = 301.68/283) with the break scale A allowed to 
vary (see Eqn. A.5). The fit functions start to diverge when ex¬ 
trapolated outside the data domain, but the data points of dif- 


^ The modified Cash Statistics (Cash 1979), Caster C-Statistics, per 
degrees of freedom. 



Fig. 9. Comparison of the convolved data points and the respective con¬ 
volving model curves for the same spectrum of Fig. 8. The orange curve 
and data points are obtained from the Maxwellian-SBPL fit, and the 
grey ones are from the COMP fit. The Maxwellian synchrotron func¬ 
tion is also overlaid (green). For display purpose, the bin size has been 
increased by a factor of 5 relative to the standard bin size. 


ferent convolving models coincide almost exactly, even when A 
is left as a free parameter. This indicates that these empirical 
functions provide good descriptions of the observed data, jus¬ 
tifying our choices of models to obtain the spectral sharpness 
angles. However, we note that the fit parameters of SBPF be¬ 
come unconstrained for a varying A, which indicates degeneracy 
in the parameter space. Therefore, we follow the catalogs and fix 
A = 0.3. 

Figure 9 repeats (in grey) the data points convolved with the 
best fit function (COMP) from Fig. 8 and shows a comparison 
with a SBPF that mimics a Maxwellian synchrotron function 
(shown in orange, overlaid on the original Maxwellian in green). 
First we separately fit a SBPF to the Maxwellian synchrotron 
model in order to obtain a curve that can be used directly in RM- 
FIT. Then we fit this Maxwellian-SBPF function to the data by 
fixing the fit parameters except for the normalization factor A. 
This demonstrates how the data points can shift under convolv¬ 
ing with a strongly differing fit function. Even though the data 
points shift, the resulting fit is significantly worse (CSTAT/dof 
= 558.00/287) than the COMP fit. This is consistent with work 
by Burgess et al. (2014), who directly convolved synchrotron 
emission spectra with photon counts and found CSTAT values 
differing typically by hundreds relative to best fit curves. 

4.2. Choice of the data domain 

A key distinction between the current work and others, is that in 
this work, we want to obtain a mathematical description of the 
peak or break curvature rather than of the whole spectrum. For 
this reason, we need to test whether our results hold up under a 
change in the data domain. The considerations when choosing 
the data domain size are (1) we want to have as many data points 
as possible, while (2) we do not want to include data too far 
away from the spectral peak, which could introduce extra cur¬ 
vature effects on the low-energy end and too many upper limits 
on the high-energy end that might pull the best fit function away 
from the data points near the peak or break, or shift the inferred 


Article number, page 8 of 16 





















































Hoi-Fung Yu et al.: The sharpness of GRB prompt emission spectra 



Fig. 10. Distributions for AG = ^catalog _ gre-flt histogram) and 

A^ieft = (dashed histogram). The values of AG and A^ieft are 

normally distributed with medians at 1.7 and 0.6 degrees, respectively. 

peak or break itself. For smaller data domains, some of the fit 
parameters can be more weakly constrained than when the full 
energy domain is used. However, as we show in Sect. 4.3, the vi¬ 
olation of the synchrotron emission model cannot be explained 
by the errors on the re-fit parameters. 

In order to find the optimal data domain size, the re-fitting 
process is repeated using different values of Fught, for spectra 
of the few brightest bursts. We find that the first-upper-limit- 
data point of the BGO detector is typically at 1.5 - 3.5 times 
the value of Ep. Therefore, in order to minimise the effect due 
to the high-energy upper limits, bright = 3.0£’p is adopted. The 
above checking process is again repeated with different values 
of Fieft- It is found that data domains smaller than (O.lFp, 3.0£’p) 
usually produce large uncertainties in the hts because the data 
are insufficient to define a definite functional shape. 

In Fig. 10, we show the differences in 9 and 0ieft between 
the catalog domain size (i.e., the full GBM energy domain from 
8 keV to 40 MeV) and the data domain size mainly used in 
this paper (i.e., O.lFp to 3.0£’p), for each spectrum. We find that 
A9 and A0ieft are normally distributed with medians at 1.7 and 
0.6 degrees, respectively. This shows that while the extra curva¬ 
ture effects contributed by the data points on the flanks lead to a 
change in smoothness, the effect is small, and limiting the data 
domain size is not strictly necessary (this is also confirmed by 
Fig. 7). 

4.3. Choice of the triangle domain 

Besides the re-fitting data domain, we also check the validity of 
the triangle domain used in the computation of 6. There are 3 
choices; (1) triangle domain > data domain, (2) triangle domain 
= data domain, and (3) triangle domain < data domain. 

Triangle domain > data domain is obviously not statistically 
sound, because we have no knowledge of how the data behave 
outside the data domain. On the other hand, we need to find a bal¬ 
ance between staying as close to the peak (or break) as possible 
and measuring a meaningful amount of curvature. As discussed 


in Sect. 4.2, the choice of Xright is already limited by the upper 
limits, so we concentrate on checking the consistency of xieft- 

In Fig. 11, we show the difference between the sharpness 
angles of the Maxwellian synchrotron function and the fitting 
models as a function of xieft (top left panel) and that between the 
Maxwellian synchrotron function and the all-BAND fit results 
(i.e., those described in Sect. 4.1 and shown in Fig. 7). The same 
plots for the differences of 0ieft are shown in the bottom panels. 
These plots are produced according to the procedure described 
below. First, a re-fit spectrum is randomly chosen. Second, we 
randomly draw new values of a and /3 from a uniform probabil¬ 
ity function characterized by the Icr errors of the spectrum, and 
a new value of 9 is computed. Third, we repeat the first two steps 
10,000 times and obtain the distributions for different values of 
xieft. We note that the plots have extrapolated below the Nal de¬ 
tector limit of 8 keV~ 0.03£’p for Ep ~ 250 keV (close to the 
median Ep time-resolved catalog value from Yu et al. in prep.), 
indicated by red vertical lines in Fig. 11. The shaded bands show 
the Icr region. 

For triangle domain choices where the lower boundary of 
the shaded band lies above 0, the difference between data and 
synchrotron theory is the clearest. The plot therefore shows how 
^left = 0.1 robustly leads to an unambiguous result. This is true 
for other choices of xieu as well, as long as xieft > 0.05. Fig. 11 
also shows that, while setting xieu = 0.3 rather than 0.1 leads to a 
larger safety margin, the difference in actual angle is negligible. 
By extrapolating the triangle domain boundary xmi to very small 
values, the long side of the triangle will eventually align with the 
left power-law asymptote. For any basic synchrotron spectrum, 
the left angle will then approach 0ieft = sin“*(3/4) a; 48.6 de¬ 
grees, corresponding to the well-known synchrotron line-of- 
death slope. Subsequently comparing this angle to one inferred 
from a best fit, therefore then becomes equivalent to testing for 
violation of the synchrotron line-of-death. We note however, that 
this analysis indicates a large error margin and extrapolating be¬ 
yond the data domain. 

5. Theoretical implications 

Our results show that for most GRB prompt emission spectra, 
an explanation in terms of synchrotron radiation can be prob¬ 
lematic. In the internal shocks of GRBs, a single-electron emis¬ 
sion function is obviously non-realistic (as there must be multi¬ 
ple electrons in the outflow) and a Maxwellian population drawn 
from a single temperature is the limiting case. Even this limiting 
case is already too wide to fit most GRB time-resolved spectra. 

The minimum variability timescale (MVT, e.g., Bhat 2013; 
Golkhou & Butler 2014) of the light curves is thought to be re¬ 
lated to the actual dynamical timescale of the emission process. 
Therefore, if the temporal bin widths of our spectra are larger 
than the MVTs, then in the time-resolved catalog we are still 
looking at averaged time-resolved spectra that are less averaged 
than those in the time-integrated catalogs. In Fig. 12, we plot 
9 against temporal bin widths per MVT (for the computational 
method of the MVT, see Bhat 2013). It is observed that, in 1,064 
spectra (49 spectra were excluded because they belong to bursts 
with no MVT due to bad or not enough GBM data), only 4.4% 
of the spectra have bin width less than the MVT for the respec¬ 
tive burst. This means that the problem for the synchrotron the¬ 
ory may be even more severe, since our spectra are smoothened 
already. However, this picture is complicated by the possibility 
that the MVT is time and energy dependent. Golkhou & Butler 
(2014) use another method to compute the MVT, which is con¬ 
sistent with our method (Golkhou et al. 2015). The uncertainties 
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Fig. 11. Top left panel: Difference between the sharpness angles of the Maxwellian synchrotron function and the re-fitting models, 6n,axw “ 
as a function of xieft. Top right panel: Same for 0maxw - ^band- Bottom left panel: Same for 6ieft,maxw - ^left,model- Bottom right panel: Same for 
^left.maxw “ ^left.BAND- The red Vertical lines show the Nal detector limit of 8 keV-^ 0.03£p for Fp 250 keV. The shaded regions show the Icr 
regions. See main text for details about the plots. 


of the MVT are worthy of independent studies which are beyond 
the scope of this paper. 

For many years, the Band function has been assumed to be 
the appropriate mathematical description in most of the GRB 
prompt spectral studies. As already shown in Sect. 4.1, the fact 
that most spectra are best fit by the Comptonized model (both 
time-integrated and time-resolved) shows that the high-energy 
tail of the prompt spectrum is actually sharper than a Band func¬ 
tion would predict (i.e., maybe somewhere in between BAND 
and COMP). In a recent study using a subsample of GBM bursts 
which occurred in the Fermi Large Area Telescope® (LAT, At¬ 
wood et al. 2009) field of view but remained undetected, Ack- 
ermann et al. (2012) showed that the Band function’s [i (as ob¬ 
tained from GBM spectral fits) is too hard to be consistent with 
the LAT upper limits. All these results are indicating that the 
Band function can lead to incorrect interpretation of the data. To 
resolve this problem, there are at least two ways: (1) to invent 
another empirical mathematical function and then again try to 


® The Fermi Large Area Telescope is a pair production telescope cov¬ 
ering the energy range from 20 MeV to 300 GeV. 


interpret the parameters of this new function by physics; or (2) 
to fit the observed spectrum directly by physical models. 

It is difficult to construct another empirical function which 
can improve upon the Band function, because it is already very 
simple in a statistical sense: it has only four parameters, and 
COMP has three. Yu et al. (2015) have shown that a triple power 
law with sharp breaks, in which the power-law indices have al¬ 
ready been constrained according to the fast- or slow-cooling 
synchrotron models, could only perform as good as the Band 
function. They have found that in many cases even an extra 
blackbody is needed to describe the spectral curvature. Recently, 
more and more studies are being performed using physical fit¬ 
ting models (e.g., Burgess et al. 2011, 2014) and simulations 
under more realistic physical conditions, e.g., varying magnetic 
fields (e.g., Uhm & Zhang 2014). However, without knowledge 
of the emission process, it is difficult to formulate a sufficiently 
well-constrained fit function. Furthermore, there may be multi¬ 
ple emission mechanisms at work, the sum of which forms the 
observed prompt spectra. 

The fitting results obtained using semi-empirical synchrotron 
models (e.g., Yu et al. 2015) and physical synchrotron models 
(e.g., Burgess et al. 2011, 2014) show that extra thermal com- 
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Fig. 12. Sharpness angles plotted against the temporal bin widths per 
MVT. Red dat points show spectra best fit by COMR blue by BAND, 
and green by SBPL. The vertical dash-dotted line shows where the bin 
width equals the MVT, only 4.4% of data points are located to the 
left of the line. The horizontal lines show the limits of the normal¬ 
ized blackbody (dotted), single-electron synchrotron (solid), and syn¬ 
chrotron emission from a Maxwellian electron distribution (dashed). 


portents are needed to fit the data. The resulting poor CSTAT 
values and systematic residual trends indicate that a pure non- 
thermal synchrotron emission function is inconsistent with the 
data at the peak or break energies, and thus cannot be the dom¬ 
inant process which contributes to the observed flux around this 
energy range. The distribution of spectral peak sharpness val¬ 
ues that we report in this paper implies that any model based on 
standard synchrotron theory without additional radiative mecha¬ 
nisms will systematically struggle to capture the spectral curva¬ 
ture of the prompt emission. This will manifest itself in relatively 
poor CSTAT values and systematic trends in the fit residuals. 

Recently, Axelsson & Borgonovo (2015) have shown that us¬ 
ing the full-width-half-maximum measurement of GRB prompt 
emission spectra taken from the BATSE 5B GRB spectral cat¬ 
alog (Goldstein et al. 2013) and 4-years Fermi GBM GRB 
time-integrated spectral catalog (Gruber et al. 2014), a signifi¬ 
cant fraction of bursts (78% for long and 85% for short GRBs) 
could not be explained by a Maxwellian population-based slow- 
cooling synchrotron function. Our results show that using the 
time-resolved spectra this violation is actually more severe, with 
over 91% of spectra obtained from long bursts violating the 
Maxwellian synchrotron function drawn from a single temper¬ 
ature, which is already a limiting case. 

As can be seen from Fig. 6, it is obvious that a small num¬ 
ber of Planck functions are not enough to reconstruct the ob¬ 
served spectral shape. From the observational point-of-view, fit¬ 
ting many blackbodies (with many parameters) is statistically 
meaningless, although maybe a sufficiently simple function de¬ 
scribing a continuum of temperatures can be formulated. On the 
theoretical side, simple photospheric models also show difficul¬ 
ties in explaining the observed data. For example, early theo¬ 
retical studies of a pure thermal origin of GRB prompt emis¬ 
sion, such as from freely expanding photospheric outflows with 
no baryonic matter or magnetic field (Goodman 1986; Paczyn- 


ski 1986), have shown difficulties in explaining the shape of 
the prompt emission phase and the two evolutionary trends of 
Ep (i.e.. hard-to-soft evolution and intensity tracking, see, e.g., 
Ford et al. 1995). Recent studies (e.g., Pe’er et al. 2006; Gian- 
nios 2008; Pe’er & Ryde 2011; Ryde et al. 2011; Vurm et al. 
2011; Lazzati et al. 2013) suggested that the Band function can 
be reconstructed from a thermal model. However, Deng & Zhang 
(2014) claim that the hard-to-soft evolution of Ep is difficult to 
reproduce under natural photospheric conditions. 

A frequently discussed alternative to the baryonic composi¬ 
tion of the jets in GRBs is a magnetically, or Poynting flux, domi¬ 
nated jet (Thompson 1994; Drenkhahn & Spruit 2002; Lyutikov 
& Blandford 2003). In this scenario, the magnetic field domi¬ 
nates the energy density in the emitting region. Thus, the domi¬ 
nant emission mechanism will be synchrotron emission from rel¬ 
ativistic electrons, since no cooling mechanism is known which 
is faster (see, e.g., Beniamini & Piran 2014). Our observational 
results therefore also pose a challenge to Poynting flux domi¬ 
nated models, although Compton up-scattering from seed pho¬ 
tons in the environment of an emerging Baryon-free jet offer a 
potential means of combining strongly magnetic outflows with 
a thermalized component or sharp spectrum (see Gill & Thomp¬ 
son 2014, for a recent example). Moreover, Beloborodov (2013) 
argues that other optically thin emission models share the same 
problems of the synchrotron emission models, e.g., pitch-angle 
synchrotron radiation (Baring & Braby 2004) when the scatter 
angle in the comoving frame is not isotropic, and jitter radiation 
in turbulent magnetic fields (Medvedev 2000). 

Finally, we compute the average time-resolved sharpness an¬ 
gles and left angles, {6} and (0ieft), weighing each spectrum 
equally. In Fig. 13, we compare (6) and (0ieft) to the sharpness 
angles and left angles computed using the time-integrated cat¬ 
alog (Gruber et al. 2014), 6'"^ and for every burst in our 
sample (listed also in Table B.l). In the left panel, green color 
indicates the 7 bursts (10%) whose average sharpness angles are 
consistent with the Maxwellian synchrotron limit (individual 9 
values can still be inconsistent, see, e.g., GRB 090829.672 from 
Fig. 5), orange color indicates the 55 bursts (79%) that are incon¬ 
sistent with the Maxwellian synchrotron limit but consistent on 
average with the single-electron synchrotron limit, and red color 
indicates the 8 bursts (11%) that are inconsistent with the single¬ 
electron synchrotron limit. Similarly, in the right panel, green 
color indicates the 13 bursts (19%) whose average sharpness an¬ 
gles are consistent with the Maxwellian synchrotron limit, or¬ 
ange color indicates the 43 bursts (61%) that are inconsistent 
with the Maxwellian synchrotron limit but consistent on average 
with the single-electron synchrotron limit, and red color indi¬ 
cates the 14 bursts (20%) that are inconsistent with the single¬ 
electron synchrotron limit. We note that the error bars of (0) and 
(0ieft) represent the standard deviations SD - yjiO^) - {9}^ and 
~ which indicate the spread of the angle 

distributions within each burst. The etTor bars of 0”' and 0|"jj are 
computed using the same procedure as described in Sect. 3, and 
are relatively small because the parameters are better constrained 
by higher photon counts. 

Figure 13 shows that the time-integrated angles are system¬ 
atically larger than the average time-resolved angles for individ¬ 
ual bursts, and the data points lie closer to (or sometimes even 
above) the diagonal in the right panel. One reason for this is that 
the spectral evolution of Ep is corrected for when computing (0), 
but not when computing 0“‘. Another reason is that rotation of 
the triangle between spectra, where a decrease in 0ieft is compen¬ 
sated for by an increase in 0right = 0 - 0ieft, or vice versa (i.e., a 
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Fig. 13. Left panel: Comparison between the average sharpness angles, (9), to the sharpness angles computed using the time-integrated catalog, 
0‘"’. Right panel: Comparison between the average left angles, (Sicft), to the left angles computed using the time-integrated catalog, The dash- 
dotted line shows x = y. The solid and dashed lines show the single-electron synchrotron and Maxwellian synchrotron limit, respectively. We note 
that the error bars of (9) and (6ieft) represent the spread in 9 and Sicft. See main text for the color-coding and details about the plots. 


specific joint change in power-law indices), is also corrected for 
when computing {0}. This latter compensation is not possible for 
(flieft), and the the data points in the right panel of Fig. 13 there¬ 
fore lie closer to the diagonal. We also emphasise that different 
light curve binning methods are used in the time-resolved and 
time-integrated spectral catalog. In our time-resolved analysis, 
as mention in Sect. 2, the light curves are binned with S /N = 30, 
and then those spectra without a peak or break are excluded. In 
the time-integrated catalog (see, e.g., Gruber et al. 2014), all time 
intervals with S/N > 3.5 are included. The fact that fewer bursts 
in the time-integrated spectral analysis are inconsistent with the 
Maxwellian limit (44 bursts, 62% for and 62 bursts, 89% for 
^left^ underlines the importance of time-resolved analysis. 

6. Summary and conclusions 

We have computed the sharpness angles 6 of the observed time- 
resolved spectra of Fermi GRBs, and compared the values to 
the sharpest cases of the synchrotron radiation theory, namely 
the single-electron synchrotron and the Maxwellian distributed 
synchrotron emission function. We find that over 91% of the 
observed spectra are sharper than the Maxwellian synchrotron 
model predicts, indicating that synchrotron radiation cannot be 
responsible for the peaks or breaks of GRB prompt emission 
spectra. No general evolutionary trend is observed for 6 within 
bursts. Moreover, the Maxwellian synchrotron function can only 
contribute up to 58)^jg% of the peak fiux. We conclude that 
the underlying prompt emission mechanism in GRBs must pro¬ 
duce spectra sharper than a Maxwellian synchrotron function but 
broader than a blackbody. 

It is still possible for synchrotron emission to dominate the 
spectrum away from the peak or break observed in the GBM en¬ 
ergy range (e.g., at the LAT energy range). Also, a sub-dominant 
synchrotron component can allow for a continuous connection 


to the afterglow phase, where synchrotron emission is typically 
dominant (see, e.g., van Eerten 2015, for a recent review). The 
transition between prompt and afterglow is then marked by the 
disappearance of the non-synchrotron (likely thermal) compo¬ 
nent. There are other theoretical possibilities to explain GRB 
prompt emission, such as the collisional model of electron- 
positron pairs (e.g., Beloborodov 2010). For recent reviews on 
GRB prompt emission mechanisms, see, e.g., Zhang (2014) and 
Pe’er(2015). 

A possibly similar inference can be made on the related phe¬ 
nomena of prompt optical emission showing a similar tempo¬ 
ral profile as the gamma-ray emission (e.g., Elliott et al. 2014; 
Greiner et al. 2014) or very early X-ray hares (e.g., Pe’er et al. 
2006; see also Hu et al. 2014 for a recent large Swift sample 
study): if the prompt emission is not dominated by synchrotron 
emission, this is likely the case for this longer wavelength emis¬ 
sion as well (see, e.g., Starling et al. 2012; Peng et al. 2014). 

We demonstrated in this paper a method to quantify the shape 
of the observed GRB spectra that provides a clear tool for dis¬ 
tinguishing between various standard emission functions. Ulti¬ 
mately, the question as to the viability of any particular emission 
model can only be fully resolved if complete spectral predictions 
for that model are tested directly against photon counts (see, e.g., 
Burgess et al. 2014). Our paper predicts that any model based on 
standard synchrotron theory without additional radiative mecha¬ 
nisms will systematically struggle to capture the spectral curva¬ 
ture of the prompt emission. This will manifest itself in relatively 
poor CSTAT values and trends in the fit residuals. 
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Appendix A: Fitting functions 

The Comptonized model (COMP) is a power-law model with a 
high-energy exponential cutoff; 


fcoMp(E) - A 



(a + 2)E 


(A.l) 


where A is the normalization factor at 100 keV in units of 
ph s“^ cm“^ keV a is the power-law index, and E^ is the peak 
energy in the vFy space in units of keV. 

The Band function (BAND) is a model which a low-energy 
cutoff power law and a high-energy power law joined together 
by a smooth transition. It is an empirical function proposed by 
Band et al. (1993): 


/band(£’) - A < 


( E Y 
( 100 keV ) 

( E -f 

V 100 kevy 


(a+2)E] 

E, 


exp 
exp (y6 - O') ( 


J- 


E <Ec 


100 keV J 


: E >E^ 


(A.2) 


where 



In Eqns. (A.2) and (A.3), A is the normalization factor at 
100 keV in units of ph s ' cm“^ keV *, a is the low-energy 
power-law index, /? is the high-energy power-law index. Zip is the 
peak energy in the vFy space in units of keV, and E^ is the char¬ 
acteristic energy where the low-energy power law with an expo¬ 
nential cutoff ends and the pure high-energy power law starts, 
in units of keV. We note that when p —> -oo the Band function 
reduces to the Comptonized model. 

The smoothly broken power law (SBPL) is a model of two 
power laws joined by a smooth transition. It was first parameter¬ 
ized by Ryde (1999) and then re-parameterized by Kaneko et al. 
(2006); 


/ F A 

/sbpl(£) = A - 10(«"«piv) (A.4) 

^ llOOkeV/ ^ 

where 

a - mAln(^^^^),apiv = mAIn^) , 

< (A.5) 

_ \og(E/Eb) _ log(lQ0keV/i-i,) 

y ‘1 — 2 ’ “pi'' ~ 2 

In Eqns. (A.4) and (A.5), A is the normalization factor at 
100 keV in units of ph s * cm“^ keV ', a and p are the low- 
and high-energy power-law indices respectively, £b is the break 
energy in units of keV, and A is the break scale. Unlike the Band 
function, the break scale is not coupled to the power-law indices, 
so SBPL is a five-parameter-model if we let A free to vary. It is 
fixed at A = 0.3 in all the Fermi GBM GRB catalogs and is 
therefore adopted in this paper. 

The peak energy of SBPL in the vFy space can be found at 

Zip = lO^Zib ,x = Atanh^* — j . (A.6) 

We note that Eqn. (A.6) is only valid for a > -2 and p < -2. 


Article number, page 14 of 16 



Hoi-Fung Yu et al.: The sharpness of GRB prompt emission spectra 


Appendix B: Comparison between the average time-resoived and time-integrated sharpness angies 


Table B.l. Comparison between the average time-resolved and the time-integrated sharpness angles. Column (1) lists the CRB names using the 
Fermi CBM trigger designation that is assigned for each new trigger detected. Column (2) lists the numbers of spectra used in averaging 0 and 
Sieft, N, for individual bursts. Columns (3) and (5) list the average time-resolved sharpness angles (0) and left angles (0ieft), and Cols. (4) and (6) 

list their respective standard deviations SD = ^j{0^) — {OY and SDieft = “ (^left)^- Columns (7) and (9) list the time-integrated sharpness 

angles and left angles and Cols. (8) and (10) list their respective errors. 


GRB name 
(1) 

N 

(2) 

{0) 

(3) 

SD 

(4) 

(5) 

(6) 

glint 

(7) 

crint 

(8) 

/Tint 

^left 

(9) 

int 

‘^left 

(10) 

080817.161 

14 

113.72 

11.40 

58.10 

6.09 

141.78 

3.11 

57.02 

0.53 

080825.593 

12 

104.97 

17.79 

50.46 

5.20 

125.47 

3.03 

50.57 

0.68 

080916.009 

12 

97.13 

7.84 

52.77 

3.85 

158.31 

8.04 

65.02 

4.35 

081125.496 

4 

79.65 

9.23 

43.98 

4.68 

125.60 

5.95 

48.77 

1.37 

081207.680 

7 

102.97 

21.57 

50.74 

7.73 

134.76 

3.85 

49.36 

0.67 

081215.784 

16 

113.42 

20.45 

51.99 

4.89 

128.58 

2.33 

53.27 

0.35 

081224.887 

11 

88.08 

11.32 

48.23 

5.65 

97.49 

0.89 

52.98 

0.44 

090131.090 

1 

125.13 

- 

66.06 

- 

149.10 

2.88 

63.36 

2.07 

090328.401 

7 

109.53 

9.37 

58.74 

4.45 

115.07 

0.90 

61.40 

0.42 

090424.592 

27 

94.86 

8.24 

51.65 

4.05 

121.41 

2.46 

59.76 

0.40 

090528.516 

9 

113.35 

8.92 

60.55 

4.21 

122.99 

1.15 

65.08 

0.53 

090530.760 

4 

74.46 

3.24 

41.34 

1.68 

122.69 

1.43 

58.23 

0.74 

090618.353 

47 

118.12 

9.18 

61.19 

4.17 

140.99 

0.93 

62.03 

0.39 

090626.189 

10 

110.62 

10.07 

55.91 

7.11 

146.27 

2.43 

63.82 

0.64 

090718.762 

7 

98.76 

7.02 

53.58 

3.43 

117.73 

1.39 

62.65 

0.65 

090719.063 

19 

80.31 

7.63 

44.33 

3.92 

96.27 

0.77 

52.38 

0.38 

090804.940 

1 

81.78 

- 

45.12 

- 

90.82 

1.23 

49.69 

0.61 

090809.978 

7 

102.46 

13.01 

51.64 

4.24 

141.06 

3.46 

54.05 

0.87 

090820.027 

82 

95.73 

10.20 

49.91 

4.23 

114.84 

1.32 

51.12 

0.23 

090829.672 

32 

130.12 

13.03 

67.80 

5.62 

153.69 

2.56 

72.83 

0.38 

090902.462 

60 

98.14 

17.93 

50.32 

6.40 

96.89 

1.05 

56.00 

0.85 

090926.181 

56 

104.08 

20.17 

50.84 

4.76 

136.77 

1.80 

55.01 

0.42 

091003.191 

11 

111.50 

21.26 

59.10 

8.73 

114.55 

1.02 

61.16 

0.48 

091010.113 

3 

87.61 

14.55 

47.94 

7.25 

114.84 

1.81 

61.30 

0.85 

091120.191 

5 

104.85 

6.26 

56.53 

3.03 

113.75 

1.39 

60.79 

0.65 

091127.976 

4 

136.45 

11.17 

71.11 

4.98 

147.43 

2.08 

63.86 

2.56 

091128.285 

5 

92.07 

6.28 

50.28 

3.11 

109.96 

1.42 

59.00 

0.67 

100322.045 

17 

101.47 

8.81 

54.31 

4.86 

150.77 

0.97 

46.10 

0.51 

100324.172 

9 

83.40 

13.16 

45.83 

6.68 

90.10 

0.83 

49.33 

0.41 

100414.097 

14 

94.27 

22.51 

46.01 

6.45 

91.59 

0.55 

50.07 

0.27 

100511.035 

2 

107.98 

3.59 

58.05 

1.71 

132.58 

1.06 

69.44 

0.48 

100612.726 

4 

98.33 

7.75 

46.38 

3.58 

120.01 

3.33 

52.36 

0.96 

100707.032 

21 

96.00 

24.53 

39.00 

4.37 

145.17 

2.17 

55.78 

0.52 

100719.989 

13 

111.32 

23.61 

45.57 

8.05 

126.63 

2.91 

52.48 

0.53 

100724.029 

40 

128.33 

21.23 

53.66 

8.31 

146.58 

1.07 

54.59 

0.22 

100728.095 

20 

87.19 

4.95 

47.85 

2.48 

111.12 

4.14 

50.33 

0.39 

100826.957 

14 

131.98 

19.54 

53.98 

4.89 

149.83 

0.91 

55.65 

0.30 

100829.876 

4 

88.51 

7.22 

48.49 

3.63 

141.16 

3.82 

52.67 

1.64 

100910.818 

5 

103.49 

10.35 

55.83 

5.01 

130.08 

4.46 

57.44 

1.01 

100918.863 

26 

101.71 

8.23 

55.00 

4.02 

99.82 

0.58 

54.12 

0.28 

101014.175 

37 

122.80 

16.43 

62.07 

6.07 

148.72 

1.09 

61.83 

0.28 

101023.951 

10 

111.86 

7.84 

58.76 

4.30 

138.19 

2.83 

59.73 

0.86 

101123.952 

30 

113.28 

18.37 

55.25 

5.31 

140.21 

1.95 

56.06 

0.29 

101126.198 

11 

116.00 

5.50 

61.82 

2.58 

129.69 

1.25 

68.13 

0.57 

101231.067 

3 

92.98 

13.39 

50.64 

6.66 

127.74 

6.73 

51.13 

1.44 

110301.214 

21 

105.55 

9.15 

54.71 

4.33 

119.12 

1.28 

56.21 

0.47 

110407.998 

4 

95.66 

9.52 

52.03 

4.68 

107.89 

1.19 

58.01 

0.57 

110428.388 

5 

78.48 

8.18 

42.80 

4.64 

101.69 

3.00 

44.70 

0.54 

110622.158 

9 

104.27 

17.11 

54.30 

3.55 

123.40 

1.64 

58.65 

0.54 

110625.881 

37 

101.15 

14.11 

51.18 

5.87 

128.70 

1.67 

53.85 

0.44 

110717.319 

12 

106.86 

11.01 

57.45 

5.26 

143.11 

3.61 

57.38 

0.69 

110721.200 

9 

123.02 

11.45 

58.32 

11.02 

160.99 

1.31 

58.56 

0.51 


Continued on next page 
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Table B.l - continued from previous page 


GRB name 
(1) 

N 

(2) 

{0} 

(3) 

SD 

(4) 

(5) 

5£>left 

(6) 

gmt 

(7) 

cr“‘ 

(8) 

/oint 

^left 

(9) 

int 

(10) 

110729.142 

6 

113.66 

10.91 

60.67 

5.08 

114.58 

1.15 

61.18 

0.54 

110731.465 

7 

90.09 

13.06 

49.21 

6.50 

113.45 

6.95 

52.88 

0.77 

110817.191 

3 

96.32 

13.85 

46.54 

3.21 

151.33 

2.06 

55.12 

1.14 

110825.102 

18 

95.10 

10.12 

51.74 

4.95 

117.89 

0.81 

62.72 

0.38 

110920.546 

11 

85.26 

18.09 

38.35 

1.06 

112.64 

2.66 

46.85 

0.38 

110921.912 

9 

109.51 

6.17 

57.30 

3.38 

128.10 

4.55 

56.72 

0.49 

111003.465 

5 

101.43 

10.22 

51.49 

4.69 

140.55 

3.51 

58.05 

0.84 

111216.389 

7 

112.41 

15.99 

57.59 

2.28 

111.91 

1.23 

59.92 

0.58 

111220.486 

9 

109.50 

8.12 

58.74 

3.86 

114.75 

0.82 

61.26 

0.38 

120119.170 

14 

107.72 

10.16 

57.87 

4.81 

129.94 

3.60 

57.56 

0.67 

120129.580 

26 

97.85 

8.05 

52.85 

4.05 

119.18 

2.69 

53.05 

0.41 

120204.054 

39 

111.53 

13.74 

59.63 

6.46 

126.38 

2.02 

63.29 

0.31 

120226.871 

9 

113.08 

11.90 

56.78 

2.98 

145.26 

2.20 

53.93 

0.56 

120328.268 

22 

125.40 

22.21 

53.35 

4.62 

145.76 

1.12 

52.85 

0.43 

120426.090 

10 

99.11 

21.79 

48.54 

1.58 

110.01 

2.42 

52.02 

0.67 

120526.303 

10 

97.10 

7.83 

52.75 

3.80 

110.52 

3.43 

51.15 

0.48 

120624.933 

18 

103.02 

16.19 

54.30 

5.48 

143.06 

4.30 

57.26 

1.26 

120707.800 

12 

130.08 

17.19 

63.49 

8.75 

145.04 

2.13 

62.19 

0.93 
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